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In this paper we seek for relevant information on the asymptotic cosmological dynamics of the 
Brans-Dicke theory of gravity for several self-interaction potentials. By means of the simplest tools 
of the dynamical systems theory, it is shown that the general relativity de Sitter solution is an 
attractor of the Jordan frame (dilatonic) Brans-Dicke theory only for the exponential potential 
U{if) oc exp:^, which corresponds to the quadratic potential R(<(i) oc in terms of the original 
Brans-Dicke field (f) = exp ip, or for potentials which asymptote to exp p. At the stable de Sitter 
critical point, as well as at the stiff-matter equilibrium configurations, the dilaton is necessarily 
massless. We find bounds on the Brans-Dicke coupling constant cubd, which are consistent with 
well-known results. 
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I. INTRODUCTION 

The Brans-Dicke (BD) theory of gravity [I| represents 
the slightest modification of general relativity (GR), by 
the addition of a new (scalar) gravitational degree of 
freedom (j), in addition to the 10 degrees which are as¬ 
sociated with the metric tensor 5^1,. Besides being the 
first physically viable modification of GR, this theory has 
been cornerstone for a better understanding of several 
other modifications of general relativity. In particular, 
the equivalence between BD theory and the f{R) theo¬ 
ries of gravity has been clearly established 0 . In contrast 
to Einstein’s GR, the BD theory is not a fully geomet¬ 
rical theory of gravity. Actually, one of the propagators 
of the gravitational field: the metric tensor, defines the 
metric properties of the spacetime, meanwhile the scalar 
field (j), which also propagates gravity, is a non-geometric 
field. The latter modifies the local strength of the grav¬ 
itational interactions through the effective gravitational 
coupling Geff oc 

While many aspects of BD theory were well-explored 
in the past (see the textbooks d, Q), other aspects were 
cleared up just recently. Thanks to the chameleon effect 
d|, for instance, it was just recently understood that the 
lower experimental bounds on the BD coupling param¬ 
eter Wbd, which were set up through experiments in the 
solar system, might not apply in the large cosmological 
scales if consider BD theory with a potential. According 
to the chameleon effect, the effective mass of the scalar 
field computed in the Einstein’s frame, depends on 
the background energy density of the environment: In 
the large cosmological scales where the background en¬ 
ergy density is very small (of the order of the critical 
density), the effective mass is also very small, so that the 


scalar field degree of freedom has impact in the cosmo¬ 
logical dynamics. Meanwhile, in the solar system, where 
the averaged energy density of the environment is huge 
compared with the one in the cosmological scale, the ef¬ 
fective mass is large, so that the Yukawa component of 
the gravitational interaction associated with the scalar 
field oc is short-ranged, leading to an effective 

screening of the scalar field degree of freedom in the solar 
system.^ 

The Brans-Dicke theory has found very interesting ap¬ 
plications specially in cosmology Q, where it has been 
explored as a possible explanation of the present stage of 
the accelerated expansion of the universe Q. The prob¬ 
lem with this is that, but for some anomalies in the power 
spectrum of the cosmic microwave background Q , at the 
present stage of the cosmic evolution, any cosmological 
model has to Mproach to the so called concordance or 
ACDM model Q. The mathematical basis for the latter 
is the GR (Einstein-Hilbert) action plus a matter action 
piece: 

-S'acdm = j “ 2A) 

+ J d‘^x\/\g\CcBM, ( 1 ) 

where £cdm is the Lagrangian density of (pressureless) 
cold dark matter (CDM). On the other hand, it has been 
known for decades, that GR can be recovered from the 
BD theory only in the limit when the BD coupling con¬ 
stant Wbd —>■ CX3. 

In the references Iffl , by means of the tools of the dy¬ 
namical systems theory, it was (apparently) shown that 
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the Jordan frame (JF) Brans-Dicke theory leads natu¬ 
ral^ to the ACDM model since, as the authors showed 
in Q , the GR-de Sitter solution is an attractor of JF-BD 
theory, independent on the choice of the self-interaction 
potential for the BD scalar field.^ The interesting thing 
is that the bounds on Wbd found in Q (wbd ~ —3/2), and 
in [13 (wbd ~ — 1), are far from the solar system-based 
experimental bound Wbd > 40000 [T5l| . 

Although the chameleon effect could (in principle) ex¬ 
plain such a discrepancy between the bounds on Wbd 
based in solar system experimentation, and those based 
in cosmological considerations, nevertheless, the bounds 
estimated on the basis of cosmological arguments: Wbd > 
120 in [I^, and 10 < Wbd < 10^ in [13 , neither are con¬ 
sistent with the ones found in the references @ and M- 
Besides, we stress that the conclusion on the existence 
of the GR-de Sitter (stable) critical point independent 
on the assumed potential in iS3, is misleading. As 
a matter of fact, in [l^ the authors seem to recognize 
that the de Sitter equilibrium configuration arises only 
for the quadratic monomial V{(j)) oc 4>^, and for the lin¬ 
eal V{(j)) (X (j) potentials. Since the estimates of are 
based on the analysis of linearized solutions which are, 
in fact, unstable perturbations around the stable GR-de 
Sitter critical point and, hence, are highly dependent on 
the choice of the initial conditions, we suspect that these 
estimates could be physically meaningless. 

In this paper we shall apply the simplest tools of the 
theory of the dynamical systems - those which are based 
on knowledge of linear algebra and of the theory of the or¬ 
dinary differential equations - to uncover the dynamics of 
cosmological models which are based in the Brans-Dicke 
theory of gravity, for several self-interaction potentials, 
in a convenient phase space. Unlike iS3 , here we shall 
explore specific potentials other than (but also includ¬ 
ing) the quadratic and the lineal monomials: V{(j)) oc 
and V{(j)) oc (j), respectively (see section ITVT) . For a better 
understanding of our analysis we shall study the vacuum 
BD theory, and the BD theory with matter, separately 
(see sections IVl and IVll respectively). In the section IVlIl 
we shall show that it is not enough that the de Sitter 
solution be a critical point of the dynamical system, in 
order for the ACDM model to be an attractor of the BD 
theory with a potential. It is a necessary condition for 
the latter, that the GR-de Sitter, i. e., the de Sitter point 
which leads to 4> = (j)Q = const., to be a stable critical 
point instead. It will be shown that, as a matter of fact, 
only (exclusively) for the quadratic monomial potential 
V{(j)) oc (p, or for potentials that asymptote to (jP, the 


^ For prior works where the de Sitter solutions are investigated 
within the frame of the scalar-tensor theories, see Ref. [a , where 
de Sitter exact and intermediate inflationary solutions are found 
for FRW models with appropriate choice of the coupling function 
^bd(</^). In a it is shown that intermediate “almost de Sitter” 
solutions might arise also when CDM is included. Other, more 
resent works on this subject, are also found 0. 


ACDM model is an attractor of the BD cosmology. 

Here we shall pay special attention to the discussion on 
the actual physical meaning - if any - of the linearized 
solutions which were used in [uni' , in order to bound 
the free parameters of the BD theory (see section IVIIII) . 
Our discussion favors the estimates based on computa¬ 
tions made at the stable equilibrium configuration, over 
the estimates based on computations performed at the 
linearized solutions. For simplicity of mathematical han¬ 
dling we shall use the dilatonic field variable instead of 
the standard BD field (j). These variables are related by 
Eq. ([3]) below. 


II. BASIC SETUP 

Here we assume the Brans-Dicke theory [I| to dictate 
the dynamics of gravity and matter. In the Jordan frame, 
which is the one considered in this paper, it is depicted 
by the following action: 


S = j d^x^\ |</i? - ^(9</))2 - 2U + 2£„|, (2) 

where = g^''V = U((/) is the scalar field 

self-interaction potential, Wbd is the BD coupling param¬ 
eter, and is the Lagrangian density of the matter 
degrees of freedom. The natural units SttG = 87r/MpL = 
c = I, are adopted. For convenience we rescale the BD 
scalar field and, consequently, the self-interaction poten¬ 
tial is also redefined: 

4> = e^,V{4>)=e^U{g,), (3) 

so that, the action m is transformed into the dilatonic 
BD action: 


S = J {R - -2U + 2e-‘"£^} . (4) 

Within the context of the (low-energy) effective string 
theory, the latter action is meant, also, to represent the 
so called “string frame” representation of the theory [T^ . 
Here we prefer, for the moment, to keep talking about 
dilatonic JF-BD theory instead of string-frame effective 
action. 

The following motion equations are obtained from (jd]) : 


Gfiu — (Wbd + I) 




V^ip + {dipY = 


2 


3 -j- 2 c(Jb 


(d^U-U) 




3 -j- 2 cjb 
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where = Rf^^- g^yR/2, and is 

the stress-energy tensor of the matter degrees of freedom. 

In this paper we shall consider Friedmann-Robertson- 
Walker (FRW) spacetimes with flat spatial sections for 
which the line-element takes the simple form: 

ds^ = —dt^ + a‘^{t)Sijdx^dx^, i,j = 1,2,3. 

We assume the matter content of the Universe in the form 
of a cosmological perfect fluid, which is characterized by 
the following state equation pm = WmPm, relating the 
barotropic pressure Pm and the energy density pm of the 
fluid, where Wm is the so called equation of state (EOS) 
parameter. Under these assumptions the cosmological 
equations ([5]) are written as it follows: 


3^2 _2,Hp + U + e-^pm, 

jy Wbd .2 , o cj • I d^pU — U 

H = +2"*’+3T^ 


2 -|- Wbd (1 -|- UI m) _ 


3 “h 2 C(Jb 


6 Pm ; 


(p -\- 3II(p p — 2 


U-d^U l-3wr, 


3 + 2a;BD 
Pm 3II “h 1) Pm — 0, 


+ 


3 “h 2 cjb 


( 6 ) 


where H = d/a is the Hubble parameter. 

Due to the complexity of the system of non-linear 
second-order differential equations ©, it is a very dif¬ 
ficult (and perhaps unsuccessful) task to find exact solu¬ 
tions. Yet, even when an analytic solution can be found 
it will not be unique but just one in a large set of them. 
This is in addition to the problem of the stability of given 
solutions. In this case the dynamical systems tools come 
to our rescue. These very simple tools give us the possi¬ 
bility to correlate such important concepts in the phase 
space like past and future attractors (also saddle equilib¬ 
rium points), limit cycles, heteroclinic orbits, etc., with 
generic behavior of the dynamical system derived from 
the set of equations m, without the need to analytically 
solve them. A very compact and basic introduction to 
the application of the dynamical systems in cosmological 
settings with scalar fields can be found in the references 

Mm- 


rules one follows when choosing these variables: (i) these 
should be dimensionless variables, and (ii) whenever pos¬ 
sible, these should be bounded. The latter requirement 
is necessary to have a bounded phase space where all of 
the existing equilibrium points are “visible”, i. e., none 
of then goes to infinity. Unfortunately it is not always 
possible to find such bounded variables. 

Besides, a certain controversy is related with the actual 
usefulness of the dynamical systems approach in cosmo¬ 
logical settings due to the apparent spurious character of 
attractors [l^, [1^. In spite of this, the dynamical sys¬ 
tems theory provides powerful tools which are commonly 
used in cosmology to extract essential information on the 
dynamical properties of a variety of cosmological models, 
in particular, those models where the scalar field plays a 
role [IQI - IM [27H3a| . In the case of scalar-tensor theo¬ 
ries, in particular the BD theory of gra vity [l|, in spite 
of several published works [9l-[lll. [3J-[4^ , the dynamical 
systems are not of common usage. In this paper, follow¬ 
ing the studies in [uni, we want to show the power of 
the simplest tools of the dynamical systems theory in or¬ 
der to extract essential information on the cosmological 
dynamics of BD-based cosmological models. 

In general, when one deals with BD cosmological mod¬ 
els it is customary to choose the following variables ^ 

El: 


_ p P'' __ '^U _ d^U 


(7) 


where the tilde means derivative with respect to the vari¬ 
able T = Ina - the number of e-foldings. As a matter of 
fact X and y in Eq. 0, are the same variables which are 
usually considered in similar dynamical systems studies 
of FRW cosmology, within the frame of Einstein’s gen¬ 
eral relativity with a scalar field matter source [l^. In 
terms of the above variables the Friedmann constraint in 
Eq. dH]) can be written as 


Notice that one might define a dimensionless potential 
energy density and an “effective kinetic” energy density 


III. DYNAMICAL SYSTEMS 

As it is for any other physical system, the possible 
states of a cosmological model may be also correlated 
with points in an equivalent state space or phase space. 
However, unlike in the classical mechanics case, where 
the phase space is spanned by the generalized coordi¬ 
nates and their conjugate momenta, in a cosmological 
context the choice of the phase space variables is not a 
trivial issue. This leads to a certain degree of uncertainty 
in the choice of an appropriate set of variables of the 
phase space. There are, however, certain - not written - 



respectively, so that the Friedmann constraint can be re¬ 
written in the following compact form: 

nf + nu + uf = 1. 

The definition for the dimensionless effective kinetic 
energy density has not the same meaning as in GR 
with a scalar field. It may be a negative quantity with¬ 
out challenging the known laws of physics. Besides, since 
there is not restriction on the sign of then, it might 
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happen that Q.u = U/ZH^ > 1. This is due to the fact 
that the dilaton field in the BD theory is not a standard 
matter field but it is a part of the gravitational field it¬ 
self. This effective (dimensionless) kinetic energy density 
vanishes whenever: 

VG -6 6 

X = - => ^ tp = - In a, 

Wbd Wbd 


or if: 


X = 0 => <p = 0 => ip = const., 

which, provided that the matter fluid is cold dark matter, 
corresponds to the GR-de Sitter universe, i. e., to the 
ACDM model. 

The following are useful equations which relate H / 
and ip/ with the phase space variables x, y and 


H 

IP 




2\/6 X — Swbd — 




3 -|- 2wbd 
2 -l- Wbd (1 -I- n 


—SVEx — 6x^ + 


3 -|- 2a;BD 

3 -|- 2a;BD 

1 — 3Wn 


- 30f, 


3 “h 2ct^B 


( 10 ) 


where is given by Eq. (jS]), and it is assumed that 
r = U/{d,pUY can be written as a function of ^ [^ : 
r = r(^). Hence, the properties of the dynamical system 
(1121) are highly dependent on the specific functional form 
of the potential U = U{p). 


IV. THE DYNAMICAL SYSTEM FOR 
DIFFERENT SELF-INTERACTION 
POTENTIALS 


In this section we shall write the dynamical system m 
for a variety of self-interaction potentials of cosmological 
interest. It is worth noticing that the only information 
on the functional form of the self-interaction potential is 
encoded in the definition of the parameter T in Eq. m- 
Hence, what we need is to write the latter parameter as a 
concrete function of the coordinate ^ for given potentials. 


A. Exponential potential 


We start with the study of the exponential potential 


U{p)=M'^ 


(13) 


Our goal will be to write the resulting system of cos¬ 
mological equations (ED, in the form of a system of 
autonomous ordinary differential equations (ODE-s) in 
terms of the variables x, y, of some phase space. We 
have: 


X 


/ 


/ 

y 


(p H 




H 

Jp ’ 


ji 

= ( 11 ) 


which, in terms of the standard BD field (j) (see Eq. GD), 
amounts to the power-law potential H(</>) = in 

the action (ED. In Eq. (USD, IVP and k are free constant 
parameters. In this - the most simple - case 


C = l- 


d^U 

U 


= 1 — k, 


is a constant, so that the system of ODE-s m reduces 
dimensionality from 3 to 2. The fact that, for the expo¬ 
nential potential T = I, is unimportant in this case since, 
as said, ^ is not a variable but a constant. 

The following plane-autonomous system of ODE-s is 
obtained: 


or, after substituting equations m into (ED, we obtain 
the following autonomous system of ODE-s: 


x' = —3a; ^1 -I- VEx — ujbdX^'^ + 




3 “h 2 cjb 




1-3r;. 


-|- [2 -|- Wbd(1 + Wm)] 


■ 30®® 


y = y 


3x ( UJbdX — 


3 -|- 2a;BD 
e + 3^ , 3y2^ 

76 


3 “h 2 cjb 


2 + CUbd (1 + Wm) 3^eff 

3 -|- 2a;BD ™ 


7 =-76a; (1 - O'(r - 1), 


(12) 


x' = —3x ^I -I- 76a; — ujbdx'^'^ 

' {x+P0) 


3 + 2u;bd 

+ [2 H- a;BD(l + u/m)] x 


1-3r;. 


3 -l- 2 cjb 


311“, (14) 


y =y 


n f 4-fc\ 3(1-0 2 

, 2 -I- Wbd (1 + Wm) ooeff 
+ 3 + 2.aBB 

where, as stated before, is given by Eq. (ED- 


,(15) 
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B. Combination of exponentials 

The combination of exponentials (M^, N'^, k and m 
are free constant parameters): 


+ (16) 
which corresponds to the BD potential 
V{(l3) = 
leads to the following 

fl -mfe- 

m = (k + m) ^ \ (17) 

where we have taken into consideration that 
, ^ l-fc + (l-m)(f)%(— 

U 1+ 

As a consequence the third autonomous ODE in the dy¬ 
namical system (EH) can be written as 


r(0 


+ {k- 1)(1 - ^)^ 


(23) 


The resulting autonomous ODE - third equation in (fT^ 
- reads 


i' = (24) 

Notice that, by setting fc = 1 and then replacing fc —>■ /i 
one recovers the ODE (|20)) for the cosh potential (fl^ . 
Working in a similar way with the sinh-like potential 


{/((p) = sinh^(/r(/j), (25) 


we obtain: 


^ = 1 — fc/r cotanh(/r(p), 


and the same 


m 


+ {k- 1)(1 - 

k(l-iY 


= — V&x \k -\-m — mk — 1 — {k + m — 2)^ — .(18) 

The particular case when , m = —k, corre¬ 

sponds to the cosh potential: 


so that the corresponding autonomous ODE is the same 
Eq. ((Ml) as for the cosh-like potential. The difference 
resides in the range of the variable For the cosh-like 
potential one has: 


U{ip) = 2M^ cosh{k(f), 
for which r(^) = fc^/(l — ^)^, and 

C' = -V6X 

C. cosh and sinh-like potentials 

The cosh-like potentials 


(19) 


( 20 ) 


1 — k^<^<l + k^ (—oo < ip < oo), (26) 

while, for the sinh-like one 

1 -I- ^ < oo, 

when —00 < </? < 0, and 

—oo < ^ < 1 — kfi, 

if 0 < :^ < oo. Here we have assumed that both k and p 
are non-negative quantities {k > 0, p, > 0). 


U{(f) = cosh^(^(^). 


( 21 ) 


where k and p are constant parameters, are also very 
interesting from the point of view of the cosmology |4lj |. 
These correspond to potentials of the following kind 


V{(/)) = M‘^(l)[cosh{ln(l)>^)f , 
in terms of the original BD field (j). We have 
d^U 


^ = 1- 


U 


= 1 — kptanh{pip), 


( 22 ) 


V. VACUUM BRANS-DICKE COSMOLOGY 

A significant simplification of the dynamical equa¬ 
tions is achieved when matter degrees of freedom are 
not considered. In this case, since = 0 =)> = 

1 -I- \/^x — Wbd x'^ , then the system of ODE-s EH simpli¬ 
fies to a plane-autonomous system of ODE-s: 


^' = -V6x(l-0"(r-l). (27) 


so that 
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In the present case one has 

=2/^ = 1 + - Wed^:^, 

= a; ^WedS^ — \/6^ f2^ + 11(7 = 1, (28) 

where we recall that the definition of the effective (di¬ 
mensionless) kinetic energy density has not the same 
meaning as in GR with scalar field matter, and it may 
be, even, a negative quantity. In this paper we consider 
non-negative self-interaction potentials U{}p) > 0, so that 
the dimensionless potential energy density 17(7 = y^, is 
restricted to be always non-negative: 17(7 = 1 + \/Qx — 
> 0. Otherwise, < 0, and the phase-plane 
would be a complex plane. Besides, we shall be interested 
in expanding cosmological solutions exclusively {H > 0), 
so that y > 0. Because of this the variable x is bounded 
to take values within the following interval: 


at first sight it appears that nothing can be said about the 
stability of this solution until the functional form of the 
self-interaction potential is specified. Notice, however, 
that since ^ = 0 at this equilibrium point, this means 
that U(i^) oc e'^, i. e., the function F is completely spec¬ 
ified: r = I. As a matter of fact, the eigenvalues of the 
linearization matrix around (0,0) are: Ai = —3, A 2 = 0. 
This means that (0,0) is a non-hyperbolic point. 

We found, also, another de Sitter solution: g = — 1 => 
H = 0, which is associated with scaling of the effective 
kinetic and potential energies of the dilaton: 


P: 

17(7 

Ai = 


1 

'\/6(I -f Wed) 



6 -I- SWb 


12 -|- 17wed + Owgj, 
4 -|- 3 wbd 


1 + Wed 


■) -^2 — 0, 


(32) 


a_<x<a+.a± = ^-(^- — -j. (29) 

This means that the phase space for the vacuum Brans- 
Dicke theory d^vac can be defined as follows: 


'I'vac = {(a:, C) : a- < x < a+) , (30) 

where the bounds on the variable ^ - if any - are set 
by the concrete form of the self-interaction potential (see 
below). 

Another useful quantity is the deceleration parameter 

q = -1- = -1 - 2\/l)X -I- SwBDa:^ 

tl^ 

^ 3(l-h\/6a::-a;BDa;^)C 
+ 3 + 2u(ee ■ ^ ^ 

Seemingly, in accordance with the results of 
without the specification of the function r(^), there are 
found four dilatonic equilibrium points Pi : (xi,^i), in 
the phase space corresponding to the dynamical system 
(I27D. The first one is the GR-de Sitter phase: 


(0,0) => x = 0 => (p = ipQ, and 

= 1 => 3iF^ = U = const., 


which corresponds to accelerated expansion q = — 1. 
Given that, the eigenvalues of the linearization matrix 
around this point depend on the concrete form of the 
function r(^), 


Ai, 2 = -- l±Wl + 


8 ( 1 -n 

3(3 + 2a;BD) 


where, as before, Ai and A 2 are the eigenvalues of the 
linearization matrix around the critical point. We call 
this as BD-de Sitter critical point to differentiate it from 
the GR-de Sitter point. 

In order to make clear what the difference is between 
both de Sitter solutions, let us note that the Fried¬ 
mann constraint (O, evaluated at the BD-de Sitter point 
above, can be written as 

e p^_3iFo + ^^^-^^3iFo-C/o, 

i. e., e~‘^pm = const. This means that the weaken¬ 
ing/strengthening of the effective gravitational coupling 
(Geff oc e~‘^) is accompanied by a compensating grow¬ 
ing/decreasing property of the energy density of matter 
Pm oc , which leads to an exponential rate o expansion 
a{t) oc e^°^. This is to be contrasted with the GR-de 

Sitter solution: diJg = Gg a{t) oc , which is 

obtained only for vacuum, pvac = Gg; pm = 0. 

The effective stiff-dilaton critical points (f7|^ = 1): 


P± : (a±, 1) ^ q± = 2 + V6a±, 

Af = 6 ^I + , A 2 = 0, (33) 

are also found, where the a± are defined in Eq. (l29ll . 

In the paragraph starting below equation (IMT) , we said 
that, seemin gly (in accordance with the results of the 
references [9l-[ll]|l. the obtained critical points are quite 
independent of the form of the function T. Notice, how¬ 
ever, that this is not true at all. For the GR-de Sitter 
point, for instance, ^ = 0, which means that 


e = i- 


d^U 


= 0 ^ C/oce‘^, 


U 
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forcing F = 1. For the remaining equilibrium points, 
^ = 1 => U = const, and F =undefined. This means 
that the equilibrium points listed above exist only for 
specific self-interaction potentials, but not for arbitrary 
potentials. Flence, contrary to the related statements in 
[3-[lll, the above results are not as general as they seem 
to be. 

Given that the critical points obtained before were 
all non-hyperbolic, resulting in a lack of information on 
the corresponding asymptotic properties, in the following 
subsections we shall focus in the exponential potential 
(fT^ : U{(p) oc exTp{k(p) ^ = 1 — k, which includes the 
particular case when 

k = I => ^ = 0 U {(f) = exp if => F = 1, 

and the cosmological constant case 

fc = 0 ^ ^=1 ^ U = M‘^, 

with the hope to get more precise information on the 
stability properties of the corresponding equilibrium con¬ 
figurations.^ These particular cases: ^ = 0, and ^ = 1, 
correspond to the four critical points obtained above. For 
completeness we shall consider also other potentials than 
the exponential. 

A. Exponential potential 

Let us investigate the vacuum FRW-BD cosmology 
driven by the exponential potential (11311 . In this case, 
since ^ = 1 — fc, is a constant, the plane-autonomous 
system of ODE-s (EH), simplifies to a single autonomous 
ODE: 


{f = const), U {f) = const., i. e., the one which occurs in 
GR and which stands at the heart of the AGDM model, 
does not arise in the general case when k =/= 1. 

Hence, only in the particular case of the exponential 
potential (IT^ with k = 1 (^ = 0), which corresponds to 
the quadratic potential in terms of the original BD vari¬ 
ables: V{(j)) = M‘^(p, the GR-de Sitter phase is a critical 
point of the dynamical system (IMl) . In this case the crit¬ 
ical points are (see Eq. dSSl)): xi = 0, a;± = a±. Worth 
noticing that a;i = 0 corresponds to the GR-de Sitter 
solution = M^exptpo, meanwhile, the x± = a±, 
correspond to the stiff-fluid (kinetic energy) dominated 
phase: = 1. While in the former case the decelera¬ 

tion parameter q = —1 — H/H^ = —1, in the latter case 
it is found to be 

q = 2 -\- V6a+ > 0. (36) 

For small (linear) perturbations e = e(T) around the 
critical points: x = Xi + e, e 1, one has that, around 
the de Sitter solution: e' = —3e => e(r) oc exp(—3r), so 
that it is an attractor solution. Meanwhile, around the 
stiff-matter solutions: 

e±(r) oc e3(2+^“±)-, 

so that, if assume non-negative Wbd > 0, the points x± 
are always past attractors (unstable equilibrium points) 
since 2 -|- \/6a- > 0. For negative Wbd < 0, these points 
are both past attractors whenever Wbd < —3/2. In this 
latter case, for —3/2 < Wbd < 0, the point is a past 
attractor, while the point a:- is a future attractor instead. 

B. Constant potential U{f) = 


{k + 2 + 2a;BD) x - ■^§(1 - k) 

I -|- 2wbd/3 

X ^I -I- \^x — lubdX^'^ ■ (34) 
The critical points of the latter dynamical system are: 




The constant potential is a particular case of the ex¬ 
ponential (IT^ . when /c = 0(^ = I=> U = const). In 
this case the autonomous ODE (IMl) simplifies: 


•\/2/3 — 2(1 -|- Wbd)^ 
3 -|- 2a;BD 


^I -I- VGx - ijJbdX^^ . (37) 


Xl 


v^(l-fc) 

+ 2 + 2u;bd 




(35) 


The critical points correspond to the following values of 
the independent variable x: 


where the a± are given by Eq. (l29l) . Notice that, since _ __, 

Xi ^ 0 (but for fc = I, in which case cci = 0 and q = v6(I -I-Wbd) 

— I), there are not critical points associated with constant 

f = fo- This means that the de Sitter phase with f = 0 Since, in this case. 


(38) 


® When the critical point under scrutiny is a non-hyperbolic point 
the linear analysis is not enough to get useful information on the 
stability of the point. In this case other tools, such as the center 
manifold theorem |42|| are to be invoked. 


H 

IP 


3 — \t^UJBDX 


3 -|- 2a;BD 
H ' 

IP 


I \/Q{1 cubd):^ 


= 0 ^ H = Ho, 


(39) 













the point Xi corresponds to BD-de Sitter expansion {q = 
— 1). At Xi the effective kinetic and potential energies of 
the dilaton scale as 

_ 6 + Swbd 

Oj/ 12 + 17 wbd + 

where, as mentioned before, the minus sign is not prob¬ 
lematic since is not the kinetic energy of an actual 
matter field. As already shown - see the paragraph start¬ 
ing below Eq. (l32l) and ending above Eq. (l33l) - this point 
does not correspond to a ACDM phase of the cosmic evo¬ 
lution, since, unlike in the GR case, in the BD theory 
the effective gravitational coupling GeS oc e~‘^ is not a 
constant and, besides, the de Sitter solution H = Hq is 
obtained in the presence of ordinary matter with energy 
density pm oc G~g. For more on this see section Ivnl 
Given that under a small perturbation (e <C 1) around 

xi: 


e(t) oc exp 


4 -|- 3cl)bd 

1 -|- LUbd ) 


this is a stable equilibrium point (future attractor) if 
the BD parameter Wbd >0. In case it were a negative 
quantity, instead, Xi were a future attractor whenever 
Wbd 4/3 and 1 c^bd ^ O- 

The critical points x± in Eq. (1551) . correspond to ki¬ 
netic energy-dominated phases, i. e., to stiff-matter so¬ 
lutions = 1, where <7 = 2-1- > 0, and, under a 

small perturbation e' = A±e, 


A± = 6 




so that, assuming non-negative Wbd > 0, the points x± 
are always unstable (source critical points). In the case 
when Wbd < 0 is a negative quantity, the point X- is 
unstable if Wbd < —4/3 (the critical point x+ is always 
unstable). 


C. Other potentials than the exponential 

The concrete form of the dynamical system (P7|) de¬ 
pends crucially on the function r(^). For a combination 
of exponentials, for instance, one has (see Eq. (IT71) h 


x' = ^-3a: -I- -I- VQx - WbdCC^) , 

= —\/&x [k m — km — 1 

-(fc-I-TO - 2) ^ . (40) 


In this case (assuming that m > k), since 




l-fc-t(l-TO)(^)%(’"-'=)^ 

1+ 


(41) 


as (fi undergoes —oo <(/3<oo=> 1 — to<^< 1 — fc. 
Hence, the phase space where to look for equilibrium 
points of the dynamical system (l40l) , is the bounded com¬ 
pact region of the phase plane (x,^), given by 

= {(x, 0 : a_ < X < <a+, 1 - to < ^ < 1 - A:} , 

where, we recall, a± = -\/3/2(l± \/l + 2a;BD/3)/a;BD (see 
Eq. (1291)). 

In the case of the cosh and sinh-like potentials, Eq. 
and (1551) respectively, one has: 

x' = ^-3x -I- 3 ^^^^^ (^1 -I- VQx - WbdX^^ , 

e' = -^x(fcV-l + 2e-?"). (42) 

The difference between the cosh and the sinh-like po¬ 
tentials is in the phase space where to look for critical 
points of (l42|l . For the cosh-like potentials one has that 
the phase space is the following bounded and compact 
region of the phase plane 

: a- < X < a+, I - kp, < ^ < I + kp) , 

while, for the sinh-like potentials the phase space is the 
unbounded region U 4)'®™*'+, where 


^va'?’ = {(x, 0 : a- < a; < a-H, I + kp < ^ < oo} , 

= {(a;,?) : a- < X < a+, -oo < ^ < 1 - fc/x} . 

A distinctive feature of the dynamical systems (HOI) 
and ()42p . is that the GR-de Sitter critical point with 
X = ^ = 0, 


PdS : (0, 0) => H = Hq, (fi = ipo, 


is shared by all of them. However, as it will be shown 
in section [VHl this does not mean that for potentials of 
the kinds p6p . (I5T]) . and (1551) . with arbitrary free param¬ 
eters, the ACDM model is an equilibrium point of the cor¬ 
responding dynamical system. As a matter of fact, only 
for those arrangements of the free parameters which allow 
that the given potential has the exponential U oc exp ip as 
an asymptote, the ACDM model is an equilibrium config¬ 
uration of the corresponding dynamical system (see the 
discussion in section Em. 


VI. BRANS-DICKE COSMOLOGY WITH 
MATTER 

In the former section we have investigated the dynam¬ 
ical properties of the vacuum Brans-Dicke cosmology in 
the phase space. Here we shall explore the case when the 
field equations are sourced by CDM, i. e., by pressure¬ 
less dust with Wm = 0, and for exponential potentials 
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-0,4 -0.? -0,2 -0.' 


-0,3 -0,2 -0.1 


FIG. 1: Phase portrait of the plane-autonomous system of ODE-s 63 for k = 1, corresponding to the exponential potential 
U{ip) oc exp(^, for different values of the BD coupling constant (from left to right): (Ubd = 0, ojbd = 1-5, oJbd = 15 and ojbd = 150, 
respectively. The GR-de Sitter point {x = 0, y = 1) is always the future attractor. The stiff-matter solution in the bottom-right 
corner in each figure is a saddle critical point in the first two cases, while it is a past attractor in the last two cases. The 
remaining stiff-matter point - the one in the bottom-left corner - is always a past attractor. 


(ITSl) only, since, in this latter case, ^ = 1 — fc, is a con¬ 
stant. This means that the relevant phase space will be 
a region of the phase plane {x,y). For this case the au¬ 
tonomous system of ODE-s (fT^ results in the following 
plane-autonomous system (see equations (fTT)) and (ITSl) ): 


x' = —3x ^1 + '/Gx — 

(.t + vTs) 


y 


3 -|- 2a;BD 
l + VGi2 + Wbd) a; 


y = y 


3x ( uj-BoX — 


"x/G (3 -f 2cU] 
A — k 




+ 


, 3(1 -fc) 2 

- y 

3 + ‘Iujbd 

2 + u;bb, 3 ff 


3 “h 


,(43) 


which has physically meaningful equilibrium configura¬ 
tions only within the phase plane: 


41mat 


{{x,y) : a- <x <a+, 

<y< 


G < y < V 1 + VGx — 


UJbdX^ 


Dgtiff . 


I 1 — x/l -b 2a;BD/3 




,0 


n'tiff : I 


/ 1 + \/l + 2wbd/3 


Of = 0; 


,01 ^ of = 0; 


\ ■x/2/3wbd 

1 


■\/6(l + Wbd) 


,0 


O®® = 


12 -b ITlObo + Glu^b 
6(1 -b Wbd)^ 


P' : i - 


\j3l^ \/A: -b 4 -b 3 wb 


k \ x/2(fc “b 1) 


O®® = 


2k^ — 3k — 8 — Gujb 


2{k + iy 


P*: - 


v^(fc-i) _ 

k 2, 2cjbd k 2 2cjb 




^ 12-6fc-6fc^ + (7-2fc-5fc^)t.BB ^ 


2(fc -b 2 -b 2 cUbd)^ 


where, in the last critical point we have defined the pa¬ 
rameter: [3 = x/l -b 2ujbd/3\/8 -b Gwbd — k{k — 2). The 
equilibrium points Pstiff and Pgtiff represent stiff-fluid so¬ 
lutions, meanwhile the remaining points represent scaling 
between the energy density of the dilaton and the CDM. 

Let us to focus into two of the above critical points: 
Pgj. and P*. As it was for vacuum BD cosmology, the de 
Sitter critical point does not arise unless k = 1. In this 
latter case (fc = 1), for the last equilibrium point in Eq. 
(H3) . one gets: 


where we have considered the facts that Of > 0 and 
y G R'^ U 0. The critical points of this dynamical system 
are: 


P, : (0,1), Q = -1 (P = Ho), of = 0, Ai .2 = -3, 

where Ai and A 2 are the eigenvalues of the linearization 
matrix around P* : (0,1). This means that, for the ex- 
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ponential potential U{}p) oc expt/?, the GR-de Sitter so¬ 
lution is an attractor of the dynamical system (H51) . This 
is evident from the FIG.[T1 where the attractor character 
of the critical point (0,1) is apparent. 

For the scaling point deceleration parameter 

is given by 

k-2 
2{k + iy 

so that, for fc = 0, which corresponds to the constant 
potential U = Uq, the BD-de Sitter solution is obtained 

q = —1 a{t) oc e~'^pm = const. 

However, since 

^ 2k^ — 3fc — 8 — Ocubd 

^ 2{k + l)^ ’ 

at fc = 0, = — (4 -I- 3 wbd), is a negative quantity, 

unless the Brans-Dicke coupling parameter falls into the 
very narrow interval —3/2 < Wbd < —4/3. Hence, for 
k = 0, but for —1.5 < ojbd < —1.33, the point does 
not actually belong in the phase space 

VII. (NON)EMERGENCE OF THE ACDM 
PHASE FROM THE BRANS-DICKE 
COSMOLOGY 

This problem has been generously discussed before in 
the reference Q . The conclusion on the emergence of the 
ACDM cosmology starting from the Brans-Dicke theory, 
seems to be supported by the existence of a de Sitter 
phase, which was claimed to be independent on the con¬ 
crete form of the self-interaction potential of the dilaton 
field in P, [I^ , and then, in Ref. [ll| the same authors 
somewhat corrected their previous claim. In this sec¬ 
tion we shall address this problem and we will clearly 
show that, in general (but for the exponential potential 
U{(p) oc e'^), the ACDM model is not an attractor of the 
FRW-BD cosmology. 

Before we go any further, we want to make clear that 
the latter statement on the non-universality of the GR- 
de Sitter equilibrium point, does not forbids the possible 
existence of exact de Sitter solutions for several choices of 
the self-interaction potential (see, for instance. Ref. H). 
What the statement means is that, in case such solutions 
are found, these would not be generic solutions, but very 
particular (unstable) solutions instead, which are unable 
to represent any sensible cosmological scenario. 

Before we start our discussion, it will be useful to state 
that a de Sitter solution arises whenever 

q = -l^H = 0^H = Ho^ a{t) oc 

This condition can be achieved even if a; ^ 0. However, 
only when 

a; = 0 ip = 0 => p = po, 


the de Sitter solution can lead to the ACDM model, 
where by ACDM model we understand the FRW cos¬ 
mology within the frame of Einstein’s GR, with a cosmo¬ 
logical constant A and cold dark matter as the sources of 
gravity. Actually, only if (^ = (/jq, is a constant, the action 
o - up to a meaningless factor of 1/2 - is transformed 
into the Einstein-Hilbert action plus a matter source: 

J d'^xy^\{R-2Uo} + 2 J 

where e‘^° = l/SnGjsf. When Lm is the Lagrangian of 
CDM, the latter action - compare with Eq. © - is the 
mathematical expression of what we call as the ACDM 
cosmological model. In the remaining part of this section 
we shall discuss on the (non) universality of the ACDM 
equilibrium point. For this purpose, in order to find use¬ 
ful clues, we shall explore first the simpler situation of 
vacuum BD cosmology and, then, the Brans-Dicke cos¬ 
mology with CDM will be explored. 

A. Vacuum FRW-BD cosmology 

In this simpler situation the de Sitter phase arises only 
if assume an exponential potential of the form 

U{p)ixexpp => V((/)=M^(/^, 

which means that ^ = 0 and F = 1, are both completely 
specified, or if ^ = 1, i. e., if 

U{p) = M^ ^ 

As a matter of fact, as shown in section|Vl for exponential 
potentials of the general form: 

U{p) = =P V{(t)) = 

with k ^ 1 and k 0, the de Sitter critical point does 
not exist. In other words, speaking in terms of the orig¬ 
inal BD variables: but for the quadratic and the lineal 
monomials, V (</) oc and V{(j)) oc (/, respectively - also 
for those potentials which asymptote to either (j)'^ or <p - 
the de Sitter solution is not an equilibrium point of the 
corresponding dynamical system. 

Now we want to show that, even when a de Sitter solu¬ 
tion is a critical point of (IWll . the existence of a de Sitter 
equilibrium point in the vacuum BD cosmology, by itself, 
does not warrant that the ACDM model is approached. 
As an illustration of this statement, let us choose the vac¬ 
uum FRW-BD cosmology driven by a constant potential 
(see subsection IV BD . In this case one of the equilibrium 
points of the dynamical system (EZl): 

Xi = 1/■\/6(l -l- Wbd) 7^ 0, 
corresponds to the de Sitter solution since 

g = —1 = 0 ^ H = Hq. 
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The tricky situation here is that, although the de Sitter 
solution {H — Hq) is a critical point of the dynamical 
system dSZl), the ACDM model is not mimicked. Actually, 
at xi, 




if = 


Vqh 

Hq 

1 + Wbd 


\/6(l + Wbd) 


= 


Hot 
1 + Wbd 


■‘fo, 


i. e., the scalar field evolves linearly with the cosmic 
time t. This point corresponds to BD theory and not 
to GR since, while in the latter the Newton’s constant 
Gtv is a true constant, in the former the effective grav¬ 
itational coupling (the one measured in Cavendish-like 
experiments) evolves with the cosmic time: 


4 + 2uj^y> ^eff Hq 

= ^-e ^ -. 

o -|- ZlUbd CefF 1 + Wbd 

Taking the Hubble time to be to = 13.817 x 10® yr (as, 
for instance, in i),i- e. , the present value of the Hubble 
constant Ho = 7.24 x 10“^^ one gets 


Geff 


--7.24 X 10"^Vr-b 

1 + Wbd 


(45) 


As a consequence of the above, if consider cosmolog¬ 
ical constraints on the variability of the gravitational 
constant [d^, for instance the ones in [dj, which uses 
WMAP-5yr data combined with SDSS power spectrum 
data: 


-1.75 X 10-12yr-i < ^ < 1.05 x lO-i^yr"!, 

Ct 

or the ones derived in Ref. [d^, where the dependence 
of the abundances of the D, ^He, '^He, and upon the 
variation of G was analyzed: 

|G/G| < 9 X 10-1®yr-i. 


from Eq. (HSl) one obtains the following bounds on the 
value of the BD coupling constant: 


Wbd > 40.37 I Wbd < —69.95, 


and 


> 79.44 I Wbd < 81.44, 

respectively. These constraints contradict the results of 
[KJI, and are more in the spirit of the estimates of 


(see, also. Ref. 0). 


point of the corresponding dynamical system. However, 
do not get confused: the above statement is not true for 
any arrangement of the free constants. 

Take, for instance, the combination of exponentials. 
The GR-de Sitter point a; = ^ = 0 entails that (see Eq. 
(ITT]) 1. either k = m = 1 ^ ^ = 0, or, for to = 1, 

arbitrary k, the point is asymptotically approached as 
tp —>■ oo if fc < 1. In the former case (A: = to = 1) the 
combination of exponentials 

U{ip) = 

coincides with the simple exponential m, uiif) = 
(M^ -I- while in the latter case (to = 1, A: ar¬ 

bitrary), assuming that A; < 1, the above potential tends 
asymptotically {ip —> oo) to the exponential U{ip) « 

For the cosh and sinh-like potentials one has (see sec¬ 
tion ITVC]): 


U{ip) = m2 ± , (46) 

where the “-I-” sign is for the cosh potential, while the 
” sign is for the sinh potential, and the 2~^ has been 
absorbed in the constant factor IVP. On the other hand, 
one has the following relationships (see section ITV C|) : 




1 — kp, 


_ g-A^V 
g/iig _|_ g-MV ’ 




1 — kpL 


gMV + g-MV 


where the left-hand equation is for the cosh-like poten¬ 
tial, while the right-hand one is for the sinh-like potential. 
Since at the GR-de Sitter point: a; = ^ = 0, then, from 
the above equations it follows that this critical point ex¬ 
ists, for the cosh and sinh-like potentials, only if kfi = 1, 
in which case, the mentioned potentials (1461) asymptoti¬ 
cally approach to the exponential as > oo: 

U{ip) fv = m 2 e‘f’. 


Summarizing: Only for the exponential potential 
U{ip) oc expy), or for any other potential which, as 
ip —> 00 , tends asymptotically to the exponential expt/?, 
the GR-de Sitter solution is an attractor of the dynam¬ 
ical system dSIl). This is easily visualized if realize that, 
by the definition of the variable 


Hence, if assume ^ = 0, which is a necessary condition for 
the existence of the GR-de Sitter point, then, necessarily: 


8 U 

—y— = 1 => U{ip) OC e'^. 


B. Other potentials 

As seen in section lV Ci for other potentials, such as the 
combination of exponentials (HU), the cosh m and sinh- 
like ()25p potentials, the GR-de Sitter solution is a critical 


C. FRW—BD cosmology with matter 

In the case when we consider a matter source for the 
Brans-Dicke equations of motion, in particular GDM, the 
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existence of a de Sitter critical point with a; = 0 => (p = 0 
- which means that the effective gravitational coupling 
is a real constant that can be made to coincide with the 
Newton’s constant - is to be associated with the ACDM 
model. 

The autonomous system of ODE-s that can be ob¬ 
tained out of the cosmological FRW-BD equations of 
motion, when these are sourced by CDM, is the one in 
Eq. (H51) . The critical points of this dynamical system 
are given in Eq. (j44]). Notice that only one of them; 

p ( v^(fc-l) P \ 

y k 2 2 ujbd k 2 2 ujbd j 

where /3 = -I- 2a;BD/3-\/8 -I- Gwbd — k{k — 2), can be 

associated with GR-de Sitter expansion, i. e., with what 
we know as the ACDM model, in the special case when 
k = 1. In this latter case P* : (0,1). Since we are consid¬ 
ering exponential potentials of the form in Eq. (US, then, 
the GR-de Sitter equilibrium configuration is associated, 
exclusively, with the potential 

^ = fc = 1 => U{tp) oc . 

Although in section I VII we have considered only ex¬ 
ponential potentials in FRW-BD cosmology with back¬ 
ground dust, it is clear that the result remains the same 
as for the vacuum case: Only for the exponential po¬ 
tential U{p) oc exp (/?, or for potentials that approach 
asymptotically to exp p, the GR-de Sitter solution is an 
equilibrium configuration of the corresponding dynami¬ 
cal system. 


VIII. DISCUSSION 

Why do our results differ from those in Ref. 
even when the tools used are the same? Although we 
should not aim here at a detailed analysis of the work 
of 1,113, nevertheless we can guess what is going on. 
To start with we shall concentrate, specifically, in the 
result related with what the authors of [3 [l3 call as 
the asymptotic value of the scalar field mass at the de 
Sitter point, which is the value of the BD scalar field mass 
computed with the help of the following known equation 

i El 113: 

= 3_^2a;BD “ ^<^^('/')] > 

or, in terms of the field variables p and U = U{ip) in Eq. 
m, the mass squared of the dilaton: 


rrP = 


3 -t- 2 ujb 


{dlU - U ), 


(47) 


According to iE3> the asymptotic value of the scalar 
field mass m\^,, at the de Sitter point, is given by 

, 1.84x10-33 

~ /q , O 

v3 -|- 2a;BD 

Then the authors constrain the BD coupling parameter 
ojbd by contrasting the above value to|* with known esti¬ 
mates signaling at to|* ~ 10“^^ eV. The obtained bound 

Wbd ~ - | + 10-^^ 

coincides with the conformal coupling value of Wbd ■ This 
value of the BD coupling parameter is a singular value 
and, if matter is taken into account, is very problem¬ 
atic since, consistency of the BD motion equations re¬ 
quire that only traceless matter can be coupled to the 
BD scalar field if Wbd = —3/2. The above bound on 
Wbd is to be contrasted with our result in the subsection 
IVII A[ which was based on the analysis of the normal¬ 
ized ratio of variation of the gravitational coupling G/G, 
which clearly excluded the possibility of Wbd = —3/2. 
Then, what is going on? 

Let us start to develop our reasoning line by recalling 
that, as properly noted in iE3, the mass squared of 
the dilaton (l47l) can be written as a function of the phase 
space variables: 


= = 


QH^y^ 
3 -|- 2cgb 


(1-C)'r(e)-1 . (48) 


If consider, for instance, the vacuum BD cosmology - see 
section IVl- one has that: 

-(l + V&x - WBDa:^) [(1 - - 1] . 

3 -I- 2a;BD ? 

Worth noticing that, since the phase space is bounded by 
the condition a- < x < a+, where the a± in Eq. (Ei) 
are the roots of the second order algebraic equation 1 -|- 
■\/6x —WbdX^ = 0, then the mass squared of the dilaton is a 
non-negative quantity, provided that r> 1/(1 — ^)^. One 
immediately notices that, at the stiff-dilaton (vacuum) 
solutions, where = 1, since = ujbbx'^ — \/6x, then 
the mass squared of the dilaton vanishes.’^ 

For the exponential potential U{(p) = exp{k(p), for 
instance, since ^ = 1 — fc, and T = 1, then 

2 6i72(p _ 1) ^ 2 

m = - 1 -I- V - ujbdX 

o + 2ujqd \ 


If introduce the dilaton’s dimensionless mass squared density 
2(fc2 i) 




3^2 3 + 2a;BD 


(^1 + VGx — UJbdX^^ , 


it is seen that it vanishes at the boundaries x = and it is 
a maximum at x = \/3/2/ujbd, where — l)/i^BD- 

The quantity is zero for the exponential potentials U± oc 

exp ± 9 ?, for which fc = dil. 


evaluated at the GR-de Sitter equilibrium point. 
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In order for the dilaton mass squared to be non-negative, 
it is required that > 1. We see that, for the particular 
case fc = 1, i. e., for the specific exponential potential 
C7((^) = M^expip, the dilaton is a massless degree of 
freedom. Hence, at the GR-de Sitter equilibrium config¬ 
uration the mass of the dilaton is necessarily vanishing. 

A general demonstration of the above statement can 
be based in Eq. gZl), where no particular considerations 
on the potential (neither on the matter content of the BD 
theory) are made. As a matter of fact Eq. (HTl) is the 
adopted definition of the mass squared of the dilaton in 
the Jordan frame [d^, which is the frame considered in 
this paper. One sees that the mass squared of the dilaton 
oc d'^U — U, vanishes provided that 

d^U = U ^ U{if) oc I U{<f) oc e^, 

where, as shown above, for the (growing) exponential 
potential U oc expy;, the GR-de Sitter critical point is 
an attractor of the corresponding dynamical system. 

Given that the scalar field is necessarily massless at 
the GR-de Sitter point, which would be the meaning of 
the tiny, yet non-vanishing, asymptotic value m|* com¬ 
puted in In this regard notice that the computations 
in are based on the linearized solutions (perturba¬ 

tions would be more precise) around the de Sitter point, 
which are valid up to linear terms in the initial conditions. 
Besides, in order to obtain the bound Wbd ~ -3/2-1-10“^^ 
on the BD coupling parameter, the authors of [1, [l^ as¬ 
sumed what they called as “special initial conditions”. 
Then, the mass of the BD scalar field computed in the 
mentioned references, is the mass of the field at the lin¬ 
earized (perturbed) solutions around the de Sitter point, 
but not at the point itself, where the dilaton is actually 
massless, as we have shown above. 

The next question would be: which is the actual mean¬ 
ing of the linearized solutions? The linearized solutions 
correspond to points in the phase space which are very 
close to the stable equilibrium point - the de Sitter criti¬ 
cal point in the present case - so that the linear approx¬ 
imation takes place: 


x{t) Ki Xc + ex(T), j/(r) Ri yc + eyir), (49) 



FIG. 2: Plots of the energy density of radiation prad = jo} 
(dash), dust pdust = la? (dots), and of the perturbed solu¬ 
tion Ppert = (dash-dots), vs the 

scale factor a, for arbitrarily chosen values of the constants: 
A = B = M = 1, ipo = Q. The plot of the energy density of 
matter in the perturbed solution is separately shown in the 
right-hand panel. 


around the GR-de Sitter point look like:^ 


x{t) k. Ae => 
^( t ) ^ 1 + B 


ip{a) Ri 


2 A 

3 


■ ‘Po, 


3(1 + 5 )^ 


These will eventually (perhaps very quickly) decay into 
the stable de Sitter solution: 


x{t) =0 => ip = ipo, 
y{0) = 1 ^ H ^ Ho = 


In the EIG. [5] the plots of the energy density of radi¬ 
ation, dust, and of the “energy density of the perturbed 
solution”, 


_ A2 _ 

Prad — T 5 Pdust — ^ 5 


Ppert — 


JVf2 

(1 + 5)^ 


where Xc , Vc are the coordinates of the given equilibrium 
point, and the perturbations e^, ^ Cy <C 1, are very small. 
These solutions can be viewed as small deformations of 
the stable GR-de Sitter solution. Just as an illustration, 
let us consider the ERW-BD theory driven by the ex¬ 
ponential potential U{p) oc exp p, in a background of 
CDM (see section Ell). The small perturbations around 
the de Sitter point P : (0,1), very quickly tend to vanish, 
restoring the system into the stable equilibrium state: 
Cx Cy oc exp(—3r), where we have taken into account 
that the eigenvalues of the linearization matrix at (0,1), 
coincide: Ai = A 2 = —3. Then, the linearized solutions 


respectively, are shown as functions of the scale factor a, 
for arbitrarily chosen values of the free constants. Notice 
that, while the energy density of “normal matter” such 
as radiation and dust, dilutes with the cosmic expansion, 
the background energy density in the linearized solution 
grows with the cosmic time, until, eventually, the per¬ 
turbed solution decays into the de Sitter expansion with 


® Here A and B are integration constants, which depend linearly 
on the initial conditions a;(0), y(0), and on other free parameters 
such as uJbo- 
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a constant (non-evolving) energy density. This means 
two things: first, that the background density of the lin¬ 
earized solution behaves as phantom matter during the 
decay into the stable equilibrium (de Sitter) state, and, 
second, that the background energy density of the lin¬ 
earized solutions is always smaller than the energy den¬ 
sity of the stable de Sitter solution. 

We can say that the linearized solutions have a small 
life-time in the sense that, within a very small amount 
of “time” T, they decay into the stable solution. It is 
then clear that the mass of the dilaton computed at lin¬ 
earized solutions, would be highly dependent on the as¬ 
sumed initial conditions, in contrast to the mass of the 
dilaton at the de Sitter point. Actually, while the mass 
of the held at perturbed (unstable) solutions, depends on 
the way the perturbations are generated, at the GR-de 
Sitter attractor - being a stable equilibrium conhgura- 
tion - the held is massless regardless of the initial con¬ 
ditions. Hence, making cosmological predictions on the 
base of linearized/perturbed solutions (around equilib¬ 
rium points) is meaningless due to the loss of predictabil¬ 
ity which is associated with the strong dependence on the 
initial conditions. The only useful information the dy¬ 
namical systems theory allow to extract from the given 
cosmological dynamical system, is encoded in the equi¬ 
librium points themselves, but not in the (linear) pertur¬ 
bations around them. The latter serve only as probes to 
test the stability of the given critical point. 

The same reasoning line applies to the computation 
of other derived quantities such as the ratio GjG (see 
section IVII A|) . In section IVlI Al we have shown that, at 
the GR-de Sitter point 


mic history into a GR-de Sitter stage, is all what one 
needs. Making definitive conclusions about the entire 
cosmic history, based in computations made at a per¬ 
turbed solution (here we are thinking in the conclusion 
on the weakening of the strength of gravity at early times 
in Ref. based in the positivity of G/G), is poten¬ 
tially misleading. The positivity of G/G at a perturbed 
solution, when compared with the negative value of G/G 
at the stable equilibrium point, may only mean that grav¬ 
ity is a bit weaker at the perturbed solution than it is at 
the stable critical point, no more. This result is closely 
connected with the fact that the background energy den¬ 
sity of the linearized solutions is always smaller than the 
energy density of the stable de Sitter solution (see the 
related discussion above). 

In a similar way we want to cast reasonable doubt on 
the conclusion in Ref. about the correspondence 

of the mean value of the BD parameter Wbd with the 
coupling parameter between the dilaton and the gravi¬ 
ton in the low-energy limit of the string effective theory 
Wbd = — 1 [Hi- As we have discussed above, in this case, 
in order to check the observational data, the authors used 
expressions for the normalized (squared) Hubble function 
{H{a)/H(ao))'^, which were computed at the linearized 
solutions (up to linear terms in the initial conditions), 
but not at the stable GR-de Sitter critical point itself. 
This means that, as previously discussed, the results of 
HO are highly dependent on the initial conditions and, 
hence, useless to make cosmological predictions. 


IX. CONCLUSION 


Geff _ Hq 
Geff 1 -|- CUbd 


—-7.24X 10-“yr-\ 

1 + Wbd 


i. e., the ratio of the variation of the gravitational cou¬ 
pling is a negative quantity (recall that we are consider¬ 
ing cosmic expansion exclusively, so that H = Hq >0). 
Hence, contrary to the result of 0, the gravitational 
coupling decreases during the cosmic expansion, result¬ 
ing in the weakening of the strength of gravity along the 
cosmic history. The explanation of the discrepancy of our 
result with the corresponding one in Ref. [l0|, is simi¬ 
lar to the explanation given above to the discrepancy in 
the bounds on Wbd- While our computation of the quan¬ 
tity G/G is done at the stable GR-de Sitter equilibrium 
point, meaning that our result is quite independent of the 
initial conditions, the corresponding computation in 0 
is done at perturbed (linearized) solutions which, as ex¬ 
plained above, are unstable and very quickly decay into 
the stable de Sitter state. The resulting computations are 
highly dependent on the initial conditions chosen to gen¬ 
erate the given perturbation around the de Sitter point. 

We hnd no reason to believe that we are living in one 
such perturbed solution and not in the equilibrium con¬ 
figuration itself (the GR-de Sitter critical point). Be¬ 
sides, if one wants to avoid the cosmic coincidence prob¬ 
lem, an equilibrium configuration which attracts the cos¬ 


In the present paper we have explored the asymptotic 
properties of FRW-BD cosmological models ([S]) , driven 
by a variety of self-interaction potentials U(ip). For this 
purpose we have used the sirnplest tools of the dynami- 
cal systems theory [9l-[lll ITgl - f^ 1 ^ . Isol - ls^ . IdOl . We 
have shown that, in spite of known results of previously 
published work @-0, tbe GR-de Sitter phase is not an 
universal attractor of the BD theory. Only for the specific 
exponential potential U{ip) oc exp:/?, which, in terms of 
the original BD field (j), amounts to the quadratic mono¬ 
mial V{(p) oc or for potentials which asymptotically 
approach to expi^ (</^)) the GR-de Sitter phase is a sta¬ 
ble critical point, i. e., a future attractor in the phase 
space. We have shown, also, that at the GR-de Sitter 
critical point, as well as at the stiff-matter equilibrium 
configurations, the effective mass of the dilaton in 
Eq. dSl), vanishes. 

We have learned that physically meaningful conclu¬ 
sions can be based only on computations performed at 
the equilibrium configurations as, for instance, at the sta¬ 
ble GR-de Sitter critical point. On the contrary, the re¬ 
sults based on computations made at perturbed solutions 
are highly dependent on the initial conditions chosen and, 
hence, useless to make physically meaningful predictions. 

In particular, the computations performed at the GR- 
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de sitter critical point yield to bounds on the value of 
the BD coupling parameter Wbd > 40.37 | Wbd < —69.95, 
or Wbd > 79.44 | Wbd < —81.44, depending on the ob¬ 
servational data assumed, which are consistent with the 
estimates of [H, [13, These results are to be con¬ 
trasted with the ones in Ref. Q: Wbd = —3/2, or in M- 
Wbd ~ —1, which were based on computations made at 
perturbed solutions. 
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